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We discuss the QCD sum-rule approach for the spacelike electromagnetic pion form 
factor in the 0(a B ) approximation. We show that the nonlocality of the condensates 
is a key point to include nonperturbative contributions to the pion form factor. We 
compare our results with the Local-Duality predictions and show that the continuum 
threshold sq(Q 2 ) parameter is highly underestimated in the Local-Duality approach at 
Q 2 > 2 GeV 2 . Using our fit for this parameter, Sq D (Q 2 ), and applying the fractional 
analytic perturbation theory, we estimate with an accuracy of the order of 1% the 0(ct 2 ) 
contribution to the pion's form factor. 
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1. A tribute to Prof. Efremov's 75th birthday 

We are happy to point out in this Festschrift, the influence of Prof. Efremov's work 
on our own research activities. 

A. V. Efremov is one of the inventors of the factorization theorems in quantum 
field theory, which form the basis of perturbative QCD applications in exclusive^ 
and inclusive^ reactions with hadrons. Without these tools, the experimental 
verification of QCD would be impossible. In cooperation with his then student 
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A. V. Radyushkin he generalized the factorization theorems for the meson form fac- 
tors, linking diagrammatic techniques with the operator product expansion (OPE). 
And all these achievements were based on previous investigations by Efremovl^tll 

Moreover, Efremov and Radyushkin have diagonalized the anomalous- 
dimensions matrix for meson operators (in leading order) in terms of Gegen- 
bauer polynomials and first obtained the asymptotic distribution amplitude (DA) 
ip(x,fi 2 — » oo) — ► ip as (x) — 6x(l — x) PJ^J Factorization theorem d^^ make it pos- 
sible to calculate various hard processes in QCD involving mesons, in which the 
meson DAs enter as the central nonperturbative input. The title of the paper^ 
" Factorization and asymptotic behaviour of pion form factor in QCD" shows ex- 
plicitly that the main interest was focused on the pion form factor (FF) for which 
the leading asymptotics at large Q 2 was revealed: F 7T (Q 2 ) — > 8 ir f 2 a s (Q 2 )/Q 2 . The 
precise value of Q 2 at which this asymptotic regime starts to prevail cannot be 
determined accurately: the estimates range from 100 GeV 2 in Refs. El [9] down to 
values around 20 GeV 2 in Refs. [10l [TTJ However, all these values are still rather far 
away from the capabilities of any operating or planned accelerator facility. 

To estimate the pion FF at intermediate Q 2 , o ne actu ally needs to e mploy some 
nonperturbative approach, like the QCD Sum Rules™ or Local Duality! * * * ' 
We present here our recent results on the pion FF obtained with QCD Sum Rules 
(SRs) with nonlocal condensates (NLCs). 



2. Pion form factor in the QCD SRs approach 

The spacelike electromagnetic pion form factor (FF) describes the scattering of 
charged particles off the pion by exchanging a photon and is defined by the following 
matrix element: 

<7T+(P')l^(0)k + (P)) = ( P+ P')^{Q 2 ) . 

Here is the electromagnetic current and q is the photon momentum q 2 = (P' — 
P) 2 = — Q 2 < in the spacelike region. To extract information about the pion form 
factor in the QCD SR approach, one needs to investigate the Axial-Axial- Vector 
(AAV) correlator of the EM current J^(x) = e u u(x)j ,1 u(x) + a d{x)^^d{x) (here 
e u = 2/3 and a = —1/3 stand for the electric charges of the u and the d quarks) 
and two axial- vector currents J§ a {x) = d{x)^^ a u{x): 

d 4 xd 4 ye^ x - p 'y^O\T[j+p(y)J»(x)J 5a (0)}\0). 

Using the standard QCD SRs t echniqu c!^^ in conjunction with the concept 
of nonlocal condensates ]^ ' l-*-8|19 | we obtain the following SR: 

sasa 

flF^Q 2 ) = J J d Sl ds 2 p 3 (si,s 2 ,Q 2 )e- (si+S2)/M2 



+ <f> G (Q 2 ,M 2 ) + <f> (qq) (Q 2 ,M 2 ), (1) 
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where M 2 is the Borel parameter, the term &g(Q 2 ,M 2 ) represents the gluon- 
condensate contribution, while the quark-condensate contribution &/q q \ (Q 2 , M 2 ) 
consists of the sum of the four-quark condensate <J>4q(Q 2 , M 2 ), the bilocal 
vector-quark condensate & 2 y(Q 2 , M 2 ), and the antiquark-gluon-quark condensate 
$> qAq {Q 2 ,M 2 ). 

The three-point spectral density is of the form 



S2,Q 2 ) 



pf\si,s 2 ,Q 2 



a s (Q 2 
An 



S2,Q 2 



(2) 



where the leading-order spectral density has been calculated long ago 



while 

the next-to-leading order (NLO) version Ap^pfai, s 2 , Q 2 ) has been derived recently 
in Ref. 1201 The contribution from higher resonances is usually taken into account 
in the form of pmi{si, s 2 ) = [1 — 6{s\ < sq)6(s 2 < sq)] Pz{si, s 2 ,Q 2 ) and contains 
the continuum threshold parameter Sq- We use in the perturbative spectral density 
the analytic version of the running coupling that avoids Landau singularities by 
construction (see for reviews in Refs. HU [55J 123) : 



a s (Q 2 ) = 



47T 



1 



A 2 

QCD 



b \ln(Q 2 /A 2 QCI} ) 
with 6 = 9 and A QCD = 300 MeV. 



A 2 

QCD 



(3) 



3. Nonlocal condensates in QCD SRs for the pion FF 



In perturbation theory, the vacuum coincides with the ground state of the free-field 
theory; hence the expectation value of the normal product is zero. Therefore, there 
are no condensate terms in perturbation theory. However, in the physical vacuum 
this is not the case. For this reason, in the standard QCD SR approach the nonzero 
quark condensate (qq) = (q~A (0)<Za(0)) appears. The value of this constant was 

19/11 

defined through comparison with experimental data for the J / tp-mesony^ Assuming 
a small coordinate dependence, the quark condensate can be represented by the first 
two terms of the Taylor expansion 



(q B (0)q A (x)) 



Sab 
4 



XAB X* 



2a s ii(qq)' 2 
81 



+ 



(4) 



where we kept the scalar and vector parts apart. Note that the condensates in this 
representation are local. 

As has been shown in Refs. [THl [T71 [551 12H the local approximation (0| is not 
reasonable for studying FFs and DAs. The reason is the unphysical behavior of 
the local condensate ^ at large x 2 , which entails a constant scalar term and a 
vector term that is even growing with the distance between the quarks x 2 . As a 
result, the nonperturbative part of the OPE linearly increases with the momentum 
Q 2 : (ci +Q 2 /M 2 ), where c\ is a dimensionless constant (not depending on Q 2 ). 
At the same time, the perturbative part decreases with Q 2 , hence generating an 
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inconsistency of the SR at intermediate and large Q 2 . Therefore, we can not rely 
upon the obtained SR for the pion FF for momentum values Q 2 > 3 GeV 2 . 

In order to improve the Q 2 dependence, one needs to modify the mo del of 
the quark-condensat e beha vior at large distances. Indeed, lattice simulations^^] 
and instanton models^EQ] indicate a decrease of the scalar quark condensate with 
increasing interquark dist ance, thu s confirming the approach of NLCsff^ The main 
strategy of the NLC SRp™ is to avoid the original Taylor expansion and 
to deal directly with the NLCs by introducing model functions that describe the 
coordinate dependence of the condensates. 

In the NLC approach the bilocal quark-antiquark condensate has the following 
formEl 

1 r°°r i 
(q A (0)q B (x)) = - / S BA (qq) f s (a) - iA x B A fv{a) e^'^da , (5) 



which, for the most general case, is parameterized by the distribution functions 
/s(a) and with = 2a s tt (qq) 2 / '81. The explicit form of these functions 

must be taken from a concrete model of the nonperturbative QCD vacuum. In the 
absence of an exact QCD solution, it was proposed^ to use the first nontrivial 
approximation which takes into account only the finite width of the spatial dis- 
tribution of the vacuum quarks: /s(a) = 8 (a — X 2 /2). This generates a Gaussian 

form of the NLC in the coordinate representation: (<lA{0)qA(x)) = (qq)e~^ A<! ^ 8 , 
which leads to the following form of the condensate contributions to the FF: 
(ci + Q 2 /M 2 ) e~ C2< 3 X q/ M , where Cj are dimensionless constants not depending 
on Q 2 . Thus, the nonlocality of the vacuum condensates generates a decreasing 
behavior of the nonperturbative part of the FF at large Q 2 . 

The same technique is applied in the case of the mixed quark-gluon condensate, 
(<7b(0)(— gA^(y) t a )qA{x)). There are two models for this condensate: the minimal 
and the improved one, see for details in Ref. 19, 26. The nonlocal gluon-condensate 
contribution produces a very complicated expression. But owing to its smallness, 
we can model the nonlocali ty of th e gluon-condensate in analogy to the quark case, 
using an exponential factor E12J5I no tably, e~ x n® ^ M . 

The described NLC QCD SR approach was used for the calculation of the pion 
FF in Ref. [26j This method yields predictions for the spacelike pion form factor (see 
Fig. [J) that compare well with the experimental data of the Cornell^ (triangles) 
and the JLab Collaborations^^ (diamonds) in the momentum region currently ac- 
cessible to experiment. These predictions cover also the range of momenta to be 
probed by the 12 GeV 2 upgraded CEBAF accelerator at the Jefferson Lab in the 
near future. This planned high-precision measurement of the pion FF at JLab will 
certainly help to check the quality of the discussed NLC models. 



a We use the Euclidean interval x 2 = —Xq — x 2 < 0. As usual in the QCD SR ap- 
proach, the Fock-Schwinger gauge is used. For this reason, all string connectors C(x, 0) = 
Vcxp[-ig 3 J*t a A1(y)dy»] = 1. 
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Fig. 1. Scaled pion form factor Q 2 F 7V (Q 2 ) for the minimal NLC model (shown as a thick broken 
line inside the shaded band delimited by the dashed lines which denote the uncertainty range). 
The improved NLC model is represented by a solid line inside the shaded band within the solid 
lines (A 2 , = 0.4 GeV 2 ). The two broken vertical lines mark the region, where the influence of the 
particular Gaussian model used to parameterize the QCD vacuum structure in the NLC QCD SRs 
is not so strong. The recent lattice result of Ref. 1311 is shown as a monopole fit with error bars 
between the two thick lines at lower Q 2 . 



4. Pion FF in the Local-Duality approach 

The LD SEPEHl is constructed from the original QCD SR in the M 2 — > oo limit. For 
this reason it has no condensate contributions. The main nonperturbative ingredient 
in this approach is the effective continuum threshold Sq D — it inherits all the 
nonperturbative information from the original QCD SR. At the (I + l)-loop order 
we have 

ss 
oo 

where 5 should be substituted by the LD effective threshold, So D '^(Q 2 ), and 
P3 \s±, S2, Q 2 ) is the three-point (/ + l)-loop spectral density. In leading order the 
integration can be done analytically and yields 



F™M(Q 2 , S) = -4t7 (l Q2 + 6S I ° 2 



4tt 2 /2 ^ Q 2 + 45 y Q 2 + 45 J 
The LD prescription for the corresponding correlatoJ^EH implies the relatic 



^ D; (°)(0)=47r 2 /'-0.7GeV 2 (7a) 

and 

C ;(1) (0) = T -^|r 7 --0.6GeV 2 , (7b) 
1 + a s (Q Q )/ir 

where Qq is of the order of Sg D '^(0). This prescription is a strict consequence of 
the Ward identity for the AAV correlator due to the vector-current conservation. In 
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principle, the Q 2 dependence of the LD parameter So D (Q 2 ) © should be determined 
from the QCD SR at Q 2 > 1 GeV 2 . But as explained in Refs. E3 GH OH the 
standard QCD SR becomes unstable at Q 2 > 3 GeV 2 because of the appearance of 
terms in the condensate contributions linearly growing with Q 2 . For this reason, this 
dependence was known only for Q 2 < 3 GeV 2 and, therefore, most authors usually 
used the constant approximation Sq D ' < "°\q 2 ) ~ Sq ''" '' (0), like in Refs. IT2| . \TE \ \TT \ 
I2TJI or a slightly Q 2 -dependent approximation Sg D '^((2 2 ) ~ in 2 f 2 /(l + a s (Q 2 )/n), 
like in Ref . U 



FF for Q 2 = 1 - 10 GeV 2 , we can estimate the effective LD thresholds Sq B (Q 2 ), 
which reproduce these predictions in the LD approach, for the two used Gaussian 
models of the QCD vacuum, the minimal and the improved one. Results are shown 
in Fig. 1. They can be represented in this Q 2 range by the following interpolation 
formulas: 



s o,min(Q 2 = x Gev2 ) = °- 57 + °- 307 tanh(0.165 x) - 0.0323 tanh(775 x) ; (8a) 



We see that Sg D (Q 2 ) in the mentioned range of Q 2 is a monotonically increasing 
function. Therefore, s^°(Q 2 ) ^ sl K (Q 2 ) rj 0.7 GeV 2 and, due to this difference, 
the LD approaches of Refs. EJ [20J [9| produce significantly lower predictions for 
Q 2 F 1T {Q 2 ) as compared with QCD SRs with NLCs. 





(Q 2 = x GeV 2 ) = 0.57 + 0.461 tanh(0.0954x) . 



(8b) 



l 



s (Q 2 ) [GeV 2 ] 



0.6 



0.8 




. 4 



0.2 



Q 2 [GeV 2 ] 



2 



4 



6 



8 



10 



Fig. 2. Effective continuum thresholds s o^ mp (Q 2 ) (solid line) and SQ^ in (Q 2 ) (dashed line) that 
approximate the NLC QCD SR results using the LD 0(a 3 (Q 2 ))-formulas. 



5. Using Fractional Analytic Perturbation Theory for two-loop 
estimates 



To estimate the Next-to-Next-to-Leading-Order (NNLO), i.e., the two-loop, con- 
tribution to the pion FF in the QCD SR approach, one needs to calculate the 
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three-loop spectral density p^\s\, s%, Q 2 ) — a complicated task. We want to avoid 
this calculation and suggest instead to use the known collinear two-loop result, the 
LD model for the soft part with an improved parameter Sg D (Q 2 ), and the match- 
ing procedure of Ref. QT| We also apply Fractional Analytic Perturbation Theory 
(FAPT) for the two-loop collinear expression in order to have an approximate inde- 
pendence with respect to the renormalization and factorization scales, see in Refs. 

HUES 

To combine the dominant, at small Q 2 < 1 GeV 2 , LD model for the soft part, 
F^ D ^°\q 2 ), with the perturbative hard-rescattering part, F%® GT> '^ 2 \q 2 ) (which 
provides the leading perturbative 0(a s ) + 0(a 2 ) corrections and is dominant at 
large Q 2 » 1 GeV 2 ), in such a way as to ensure the validity of the Ward identity 
(WI) F^ 1 '^^) = 1, we apply the matching procedure, introduced in Ref. [Til 

F WI;(2) (g2) = F LD,(0) (g2) + ( F PQCD,( 2 ) (Q2) (g) 

with Sg 2 -* ~ 0.6 GeV 2 . To test the quality of the matching prescription given by E> 



©, we compare it with the LD model ([6]) evaluated at the 0(a s )-approximation!2Q12l 
To this end, we construct the analogous 0(a s )-model F^" ;(1) (Q 2 ), where we substi- 
tute FP QCD ' (2) (Q 2 ) by F£ QCD,(1) (Q 2 ) = 2 a s (Q 2 ) s^ m {0)/n Q 2 and employ the 
same prescription for the effective LD threshold as in Refs. |H1 i.e., Eq. (|Tb[) . The 
key feature of this matching recipe is that it uses the information on F^iQ 2 ) in the 
two asymptotic regions: 

(1) Q 2 — ► 0, where the Ward identity dictates F n (0) = 1 and, hence, F^{Q 2 ) ~ 
F. LD < (0) (Q 2 ), 

(2) Q 2 -> oo, where F^(Q 2 ) ~ F% QCIU1) (Q 2 ) 

in order to combine properly the hard tail of the pion FF with its soft part. Numeri- 
cal analysis shows that the applied prescription yields a pretty accurate result, with 
a relative error varying in the range 5% at Q 2 = 1 GeV 2 to 9% at Q 2 = 3 — 30 GeV 2 . 
Note here that this recipe was proposed without the knowledge of the exact two-loop 
spectral density — this appeared laterPHl 

Now, knowing the exact expression for the spectral density (si, S2, Q 2 ), we 
can improve the representation of the LD part by taking into account the leading 
0(a s ) correction in the electromagnetic vertex. To achieve this goal, we suggest the 
following improved WI model: 

f™£\q 2 ,s) = f^\q>,s)+ s asiQ2) ( 23 



T: """ 1 ! ' " 4tt 2 / 2 tt \2S + Q 

S 



(^) do) 



with the subsequent substitution S — > Sg D (Q 2 ). Numerical evaluation of this new 
WI model in comparison with the exact LD result in the one-loop approximation 
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shows that the quality of the matching condition is improved: the relative error is 
reduced, reaching only 4% at Q 2 = 1 — 10 GeV 2 . 



We construct the two-loop WI model F^ 1 (Q 2 



,Sq D;(2) (Q 2 )) for the pion FF 



to obtain 

F^(Q 2 ,S) = F^(Q 2 ,S)- 



S a s (Q 2 ) 



4tt 2 / 2 
S 



25 



2S + 

F FAPT,(2) (Q2) 



Q 2 



(11) 



4tt 2 / 2 * y ^ ' \2S + Q 2 , 

where F PAPT '^ 2 \q 2 ) is the analyticized expression generated from F%'~ lCD '^ (Q 2 ) 
using FAPT (see Refs. EH [38l E2]) to get a result which appears to be very close 
to the outcome of the default scale setting = //| = Q 2 ), investigated in detail 
in Ref. [UJ in the APT approach. FAPT is needed here in order to obtain analytic 
expressions for the pion FF, using two possible options for the factorization scale: 

(i) For /i| = Q 2 , there appear factors of the type [a s (Q 2 )] 1 ' with fractional powers 
v = 7„/ (2 bo) due to the evolution of the pion distribution amplitude; 

(ii) For /ip = const., the factor [a s (Q 2 )] ln(Q 2 /^|) appears. 



In any case, the NNLO correction involves the analytic image of the second power 
of the coupling, A2(Q 2 )- For this reason, we call the whole F^ APT '^ 2 \q 2 ) term the 
0{Az) contribution. 



Q 2 F W (Q 2 ) [GeV 2 ] 
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Fig. 3. We show as a narrow dashed-dotted strip the predictions for the pion FF, obtained in 
the two-loop WI model, Eq. (TIT}, using the minimal (left panel) and the improved (right panel) 
Gaussian models. The width of the strip is due to the variation of the Gegenbauer coefficients a,2 
and a4 (needed to calculate the collinear part pPQ CD '( 2 ) (Q 2 )) in the corresponding shaded bands 
for the pion DA (indicated by the central solid line). Note that this dashed-dotted strip shows the 
effect of the 0(^2) correction only for the central solid curve of the shaded band. 



It is interesting to note here, that in the case of the one-loop approximation, the 
relative error of the WI model (fT0|l appears to be of the order of 10%. The relative 
weight of the 0(a 2 )-contribution to the pion FF is of the order of 10%, as has been 
shown in Refs. [UJ [22l Hence, the relative error of our estimate is of the order of 
1% — provided we take into account the O (a .^-correction exactly via the specific 
choice of s Q (Q 2 ), as done in Eq. (|8a| . 
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The results obtained for the pion FF with our two- loop model, i.e., Eq. (Ilip . 
and using the effective LD thresholds Sq D (Q 2 ), are displayed in Fig. 2. We see from 
this figure that the main effect of the NNLO correction peaks at Q 2 > 4 GeV 2 , 
reaching the level of 3 — 10%. 

6. Conclusions 

We presented here the results for the spacelike pion form factor obtained within 
the QCD SR approach, using two different Gaussian NLC models. These NLCs 
entail the decay of the nonperturbative OPE terms at large Q 2 . These NLCs make 
the QCD SR stable and enlarge the region of its applicability towards momenta as 
high as 10 GeV 2 . The principal ingredients of our approach are, besides the NLCs, 
the 0(a s ) spectral density, and the analytic Shirkov-Solovtsov coupling which is 
free of Landau singularities. Our predictions for the pion FF in the momentum 
range up to 10 GeV 2 are in a good agreement with the existing experimental data 
of the Cornell^ and the JLarj22l collaborations, as well as with a recent lattice 
calculation.^!] 

We also showed here that the LD model for the pion FF suffers from the threshold 
sq(Q 2 ) uncertainty. We fixed this uncertainty by dema ndin g that the LD model 
should reproduce the results of the Borel SRs with NLCs!^ 21 Our results show that 
So D {Q 2 ) grows with Q 2 . 

We also proved that the rough model for the matching function in Ref. IT] appears 
to be of a rather good quality (« 10%). We improved it here to reach the quality of 
« 5%. Using FAPT, and the improved matching function, wc estimated the NNLO 
correction to the pion FF to be of the order of ~ 3 — 10%. 
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